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ABSTRACT 


The Parkes Half-Jansky Flat Spectrum Sample contains a large number of 
sources with unusually red optical-to-near-infrared continua. If this is to be 
interpreted as extinction by dust in the line-of-sight, then associated material 
might also give rise to absorption in the soft X-ray regime. This hypothesis is 
tested using broadband (0.1-2.4 keV) data from the ROSAT All-Sky Survey 
provided by Siebert et ah (1998). Signihcant (> 3a conhdence level) correlations 
between optical (and near-infrared)-to-soft X-ray continuum slope and optical 
extinction are found in the data, consistent with absorption by material with 
metallicity and a range in gas-to-dust ratio as observed in the local ISM. Under 
this simple model, the soft X-rays are absorbed at a level consistent with 
the range of extinctions (0 < Ay < 6 magnitudes) implied by the observed 
optical reddening. Excess X-ray absorption by warm (ionised) gas, (ie. a ‘warm 
absorber’) is not required by the data. 

Subject headings: dust,extinction — quasars: general — X-rays: ISM — ISM: 
general 
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1. Introduction 

There have been numerous studies reporting the presence of soft 2keV) X-ray 
absorption in excess of that expected from the galaxy towards radio quasars. However, such 
studies have found very little evidence for associated optical reddening by dust (Elvis et 
ah 1994). Strong evidence for associated Mg II and soft X-ray absorption in a number of 
radio loud quasars has been conhrmed (eg. Mathur et ah 1994 and Mathur 1994), though 
evidence for associated optical extinction in these sources is weak. The lack of associated 
optical extinction may be due to a selection bias. Since X-ray absorption estimates are 
derived from spectra which require relatively large X-ray counts, such studies may be biased 
against those sources with low counts due to strong X-ray (and hence optical) absorption. 
It is also possible that X-ray absorption in most sources is dominated by the presence of 
warm (ionised) absorbers close to the primary continuum source (Pan, Stewart & Pounds 
1990). These have little broadband effect in the optical bandpass but may strongly absorb 
X-ray emission. 

If strong optical extinction by dust is known a-priori however, then associated 
absorption of soft X-rays is expected to be present at some level. Evidence for excess soft 
X-ray absorption in a number of optically reddened radio quasars has been presented by 
Kollgaard et ah (1995). These authors claim that their results are strongly model dependent, 
and are also consistent with explanations other than absorption by associated gas and 
dust. Their sample size is also too small to draw any reasonable conclusions. Evidence for 
an association was recently reported by Puchnarewicz et al. (1996) for a large sample of 
Seyfert Is and quasars. These authors found a correlation between optical spectral slope 
and optical-to-soft X-ray continuum slope that was consistent with absorption by dusty 
“cold” gas with column densities > lO^^cm”^ and approximately Galactic dust abundance. 

In this paper we consider the Drinkwater et al. (1997) sample of flat-spectrum radio 
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sources (hereafter Parkes quasars), a subset of which are known to be optically reddened 
( [Webster et ah 1995| , [Francis et ah 199^) . Unfortunately, very little spectroscopic X-ray 
data exists for these quasars to constrain absorption gas column densities. Siebert et ah 
(1998) were able to measure power-law photon indices for 105 sources in the sample. They 
did not hnd any evidence for excess soft X-ray absorption in the red sources, but they did 
hnd that (at low redshifts) the red quasars had lower soft X-ray luminosities than the blue 
quasars and concluded that the X-ray data support the importance of dust in the reddest 
quasars. 


We re-analyse the soft X-ray (0.1-2.4keV) ROSAT broadband measurements from 
Siebert et ah (1998) using a different approach. We consider a subsample of 119 out of the 
323 sources in the Drinkwater et ah sample, all of which have contemporaneous Bj and 
Kn photometry, useful for estimates of spectral slopes. About half of these sources are 
detected in soft X-rays and upper limits are available for the remainder. Our analysis uses 
a simple gas-dust absorption model to explore various correlations involving optical (and 
near-IR)Ao-soft X-ray continuum slope and optical reddening. We then use the broadband 
soft X-ray data to search for such correlations and determine whether soft X-rays are 
absorbed at a level consistent with the observed optical reddening. 


2. Quasar Sample and X-ray Data 

The Parkes Half-Jansky Flat Spectrum Sample contains 323 sources and is described 
in detail by Drinkwater et ah (1997). The sample was initially selected from the Parkes 
2.7 GHz Survey (Bolton, Savage & Wright 1979 and references therein) according to 
the following criteria: 2.7 GHz radio flux densities / 2 . 7 GHz > 0.5 Jy, radio spectral 
indices a^j < 0.5 (where fu oc z/“"). Galactic latitude |6| > 20° and B1950 declinations: 
—45° < 6 < 4-10°. By selecting flat-spectrum radio sources at high frequency, one is biased 
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towards core-dominated quasars, since lobe dominated quasars and radio galaxies have 
steeper radio spectra and hence are likely to comprise the majority of detections in low 
frequency surveys. On the basis of spectroscopic identihcation alone, the present sample 
contains a much higher quasar fraction (;^ 85%) than any existing radio sample, with a 
broad and flat distribution in redshift to z ^ 4. 

For the purposes of this paper we have used sources with contemporaneous photometric 
measurements in Bj,V, R, I, J, H, Kn, which are complete for a subsample of 119 ‘non- 
extended’ sources. These were obtained by P. Francis (Private communication; see Whiting 
et al. 1998) in April, July and September 1997 using the ANU 40" and 2.3m telescopes. 
These measurements provide reliable estimates of broadband colours, likely to be unbiased 
with respect to uncertainties from intrinsic variability. RMS uncertainties based on noise 
statistics in these magnitudes are < 0.2 mag. 

Since we are primarily interested in the properties of quasars, only spatially unresolved 
sources in the optical and near-infrared have been considered. All sources appearing 
extended in Bj and Kn, or with optical spectra showing features characteristic of those 
seen in normal nearby galaxies are excluded. We also ensured that the sources have broad 
emission lines in their optical spectra with velocity widths > 2000km s“^ (at FWHM), 
typical of normal QSOs. Redshifts are available for all the sources and span the range 
0 < z < 3.9, most of which are from new spectroscopic observations (see Drinkwater et al. 
1997). 

X-ray data for all Parkes quasars has been provided by Siebert et al. (1998). The fluxes 
are in the soft X-ray band 0.1-2.4 keV, and most were determined from the ROSAT All-Sky 
Survey and (for 49 sources) from archival pointed PSPC observations. Of the 323 sources 
in the Drinkwater et al. sample, 163 were detected in soft X-rays at the 3(T level. For 
the remaining 160 sources, 2a upper limits to the counts were determined. We have used 
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the total broadband fluxes as computed by Siebert et al. from the count rates, corrected 
for Galactic absorption only. Where available, we also use their estimates of the photon 
energy indices F (where fi, oc determined from explicit power-law fits to spectral data 

or from hardness ratio techniques (see § |4.1|) . In total, for our subsample of 119 Parkes 
quasars with contemporaneous photometry, 57 are detected in soft X-rays and upper limits 
are known for the remaining 62. 


3. Optical Extinction and X-ray Absorption Modelling 

As claimed by Webster et al. (1995), if the large spread in optical-to-near-IR colours 
of Parkes quasars is due to reddening by dust, then it is expected that the reddest quasars 
may also be absorbed in soft X-rays. We test this hypothesis by making simple predictions 
involving optical-to-soft X-ray (osx) and near-IR-to-soft X-ray {aKx) continuum slopes. 
This section will briefly outline our assumptions and predictions using a simple gas-dust 
absorption model. 

The degree of X-ray absorption by metal enriched gas primarily depends on the total 
column density of gas in the line-of-sight. In the case of the galactic ISM where the metal 
abundance is typically ^ 1% relative to hydrogen by mass (Grevesse & Anders 1991), 
hydrogen and helium are responsible for almost all of the absorption at energies ^ 2keV. We 
predict the amount of X-ray absorption expected for a given optical dust extinction measure 
by assuming for simplicity, the range in Galactic gas-to-dust ratios derived empirically from 
Ly-a absorption measurements in the Galaxy (Bohlin et al. 1978): 

A(HI + R 2 )tot ^ (5.8 ± 2.5) x lO^^ ( cm’^, (1) 

V mag J 

where Eb-v is the extinction (colour excess) m. B — V colour. This relation is also consistent 
with empirical estimates of the dust-to-gas ratio in the SMG and LMG by Bouchet et al. 
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(1985) and Fitzpatrick (1985). 

The effective X-ray optical depth, rx, defined snch that the change in finx dne to 
absorption is exp(— tx), is given by: 


rx = aEN{m + H2 ) tot ) 


( 2 ) 


where represents the effective absorption cross-section per H atom at energy E and 
iV(HI -I- H 2 )toi is defined by eqnation (|lD. We adopt cross-sections for X-ray absorption as 
derived by Morrison & McCammon (1983) for a gas with galactic ISM metal abnndances. 
Typically, cxikev — 2.42 x 10“^^cm^ and cT 2 keV — 4.30 x 10“^^cm^. Together with the ratio of 
total to selective extinction Ay/E b-v — 3.05 ± 0.15 as given by Whittet (1992), the optical 
depths at IkeV and 2keV from eqnation (^ can be written respectively: 



( 3 ) 



TakeV ^ (0.08 ± 0.04) 


As a simple estimate, if a sonrce at redshift z = 1, (the median redshift of the Parkes 


qnasar sample) snffers an intrinsic extinction Ay = 2 mag (a typical mean valne for Parkes 


qnasars; Masci 1997), then absorption by associated nentral gas wonld rednce the observed 
IkeV finx by abont a factor 1/exp(—r(i+^)ikev) ~ 1-2. For comparison, assnming the 
generic 1/A dnst law, an intrinsic extinction of Ay = 2 mag will resnlt in a decrease of the 
observed optical finx by a factor > 50. Thus, we see that the decrease of observed optical 
flux by dust extinction is generally larger than the corresponding decrease in soft X-ray flux 
by associated gas. Even for low 2 ; sources, the discrepancy in flux reductions in these two 
regimes is about a factor of 5, and increases rapidly with 2 ; due to the energy dependence 
of aE- This property will become important when we examine the correlation involving 
optical-to-X-ray continuum slopes in § |3T] . 




3.1. asx axx versus Optical Extinction 


Based on optical and X-ray spectral properties of Parkes quasars, both Webster et al. 
(1995) and Masci (1997) claimed that in most cases, the reddening is likely to be due to 
dust intrinsic to the quasars. Thus for simplicity, the following analysis will assume all 
absorbing material to be located at the redshift of the quasar. We also assume that all 
predicted quantities discussed below (ie. slopes and extinction measures) refer to the quasar 
rest frame. 


If a source with some intrinsic (unabsorbed) optical-to-X-ray power-law slope, say 
C(BXi (where oc z/“") suffers intrinsic absorption by dusty gas, then in its rest frame, the 
resulting slope (aBXo) can be written in terms of extinction optical depths as follows: 


OiBXo 


OlBXi + 


'Of — tb 

\n{vx/vB)' 


( 4 ) 


We shall also consider the near-infrared-to-X-ray continuum slope (oxXo), so that 
analogously. 


OiRXo = OiKXi + 


"Of — tk 


( 5 ) 


In {ux/i^k) ' 

According to the available data, our anlysis will assume the -bandpass (A = 475nm) for 
the optical and the iP„-bandpass (A = 2.15/rm) for the near-infrared. We shall consider 
the X-ray flux at IkeV, so that the optical depth tx in equations (^) and (^ is dehned 
in terms of the optical extinction Ay by equation (Q). Using the extinction coefficients 
Rx = Ax/Eb-v from Savage & Mathis (1979) in the optical and Whittet (1988) in the 
near-infrared for diffuse galactic dust, we hnd 


Tb ~ (1.223 ± 0.015)Ay and tk ^ (0.084 ± 0.030)71^. (6) 

By combining the above relations, we can see from equations (§) and (^) that for given 
intrinsic (unabsorbed) values aBXi and axXi, absorption by dusty-gas will predict a specihc 
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correlation between optical extinction (Ay) and the corresponding absorbed continuum 
slopes. 

To make some illustrative predictions applicable to radio quasars, the intrinsic slopes 
ttBXi and ttKXi need to be specihed. As a working measure, we assume for simplicity, 
these slopes to be those found for optically-selected quasars. As argued by Masci (1997) 
and Francis et ah (1998), this choice is based on the claim that optically-selected quasars 
are expected to be strongly biased against signihcant absorption by dust. This indeed is 
consistent by a number of studies which hnd a relatively small scatter in the optical-to- 
X-ray (IkeV) flux ratio of optically-selected quasars (eg. Kriss & Canizares 1985, Wilkes 
et ah 1994, La Franca et ah 1995). The distributions in this flux ratio indicate a mean 
value {aBXi) ~ 1-3 with dispersion a ~ 0.2, and which does not signihcantly differ between 
radio-loud and radio-quiet quasars (Wilkes & Elvis 1987, Green et ah 1995). This value 
is also consistent with that implied by a composite SED for radio loud (optically-selected) 
quasars derived by Elvis et ah (1994). Due to the absence of sufficient near-infrared data 
for optically selected samples, we adopt the iF-to-lkeV continuum slope as indicated by 
this composite where {cxkx^ ~ 1, and assume a scatter similar to that found for aBXi in 
the studies above. From these studies, our predictions will assume the following ranges: 

1.1 < aBXi < 1 - 6 ) ( 7 ) 

0.8 < axXi < 1-3. 


The radio-selected Parkes quasars for comparison show a dispersion in aBx almost hve 
times greater (see § than those of optically-selected quasars (eq. 0). Siebert et ah 
(1998) have shown that this is consistent with dust-gas absorption, however in addition, 
the larger scatter could also include effects from any of the following: an additional X-ray 
emission component over and above that of radio-quiet quasars with a wide distribution 
of strengths and/or slopes, enhanced X-ray variability, or a strong angle dependence of 
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the observed X-ray emission due to beaming. Although such effects can be important in 
some sources and contribute to the scatter, our analysis here is purely concerned with the 
hypothesis that it is mostly due to variable amounts of absorption in the optical and X-ray 
bands. 

Given the ranges in intrinsic slopes defined by equation (|^) and the ranges defined by 
the uncertainties in the extinction measures rikev (eq. [@) and tb, tk (eq. p), we show in 
Figure 0 the expected range spanned by ubXo axxo a function of Ay- There appears 
to be a distinct behaviour in each of these slopes with Ay, with an anti-correlation in aBXo 
and correlation in axXo- As discussed in § |^ and as seen in equations (^) and (]^), this is 
due to the amount of absorption expected in Bj and Kn relative to that at IkeV. From the 
above discussion, we have typically: tb/tx ~ 2.7 and tkItx ~ 0.2, so that the correlation 
involving axxo is purely due to the fact that absorption causes a greater decrease in soft 
X-ray flux than that in the near-infrared. It is important to note that the trends of these 
correlations sensitively depend on the assumed dust-to-gas ratio (eq. [|l|) and metallicity 
(ie. the cross-section ge in eq. [0), which for simplicity were fixed to local ISM values. A 
variation in either of these quantities by at least a factor of five is likely to change the sign 
of the correlation. 

The predictions in Figure |l| are represented by the following linear relations: 

= (1.3;°i) - (0.13;°:°3) X Ay (8) 

+ (0.05;°:°^) X Ay. 

Given the uncertainties, the slope of the olkx versus Ay relation is insignificant, however 
for typical galactic ISM conditions, the crucial feature in Figure |l| is the opposite trend 
predicted in each of these slopes in the presence of absorption by dusty gas. This will 
provide a powerful diagnostic for testing the dust reddening hypothesis for Parkes quasars. 
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EDITOR: PLACE FIGURE 0 HERE. 

4. Comparison with Data 

This section will compare the predictions of Fignre with estimates of the slopes 
involving the soft (0.1-2.4 keV) X-ray bandpass {obx and axx) derived from the available 
ROSAT data, and extinctions {Ay) derived from the observed optical reddening. To 
facilitate a direct comparison with the predictions of Fignre |l| however, we first transform 
these qnantities into the sonrce rest frame. 


4.1. Spectral Slopes and Transforming to the Rest Frame 

Dne to the freqnency dependence of optical extinction and X-ray absorption, a similar 
dependence of these qnantities on the redshift of the absorbing material in an observer’s 
frame is expected. This implies that a plot similar to Fignre |I] which uses observed quantities 
will also include a hidden and complicated dependence of the absorption on redshift. The 
added effects of changing spectral shape with source redshift (ie. K-correction effects) 
will also introduce ambiguities. Thus, assuming the absorbing material is intrinsic to the 
quasars, we will transform all observed quantities to the source rest frame. This will enable 
a direct and unambiguous comparison with the predictions of Figure |I]. 

We have used the total 0.1-2.4 keV X-ray fluxes as computed from the ROSAT count 
rates by Siebert et ah (1998), which were corrected for Galactic absorption only. The 
observed spectral indices (ie. aBx(obs) oc log [/a,//^]), were computed by first deriving 
monochromatic X-ray fluxes fx at IkeV. These were determined from the total broadband 
fluxes and available upper limits assuming a power-law continuum oc between 
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0.1-2.4 keV. Where available, we have used the photon indices F (where ax = F — 1) 
determined from explicit power-law hts to spectral data and hardness ratio techniques by 
Siebert et ah (1998). These were available for 71 of the 119 sources. When an individual 
photon index was not available, the average photon index for radio-loud quasars, ax = 1 
(eg. Schartel et ah 1996) was used. 

The rest frame spectral indices, for example aBx(rest), were determined by applying a 
simple K-correction which assumes a power-law in each of the optical and soft X-ray bands. 
If as and ax are respectively the optical and X-ray power-law slope, z the source redshift 
and asx(obs) the observed optical-to-X-ray slope, then the corresponding value in the rest 
frame is given by: 

{aB - ax)\og{l + z) . I .. 

aBx(rest) = -^^- - —r-h asx(obs). (9) 

log {yx/vB) 

X-ray slopes ax are taken from Siebert et ah (1998) as discussed above, and optical slopes 
as for each source were determined from the contemporaneous photometry, measured 
between the Bj and / (~ 0.9/im) passbands. In our determination of the near-IR-to-X-ray 
slope as'x(i'est), the required near-IR slopes ax were determined between the and 
passbands. 

Rest frame optical extinctions Ay in each source were derived from the observed 
optical-to-near-IR reddening as dehned by the contemporaneous colours Bj — I and Bj — 
and assuming intrinsic (unabsorbed) colours as measured in optically-selected quasars. As 
discussed in § ^Tj optically-selected quasars are expected to be strongly biased against 
signihcant absorption by dust. This claim is consistent with their relatively small scatter 
in colours (eg. Francis 1996), which he predominately on the blue tail of the Parkes-quasar 
colour distribution. From quasi-simultaneous optical/near-IR photometry by Francis (1996) 
of a subset of 37 quasars drawn from the optically-selected LBQS sample of Hewett et al. 
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(1995), we find the following mean values and dispersions in intrinsic colours: 


{B - K)i = 2.3 ± 0.5 and {B - I)i = 0.9 ± 0.4. (10) 


We assume these values represent the intrinsic (unabsorbed) colours of Parkes-quasars. 
For our redshift range of interest: 0 < z < 3, we also find that for these optically-selected 
quasars, the colours show no signihcant dependence on redshift. We therefore assume the 
intrinsic colours (eq. [|^) to be independent of redshift in an observer’s frame. 


Given a general dust extinction curve dehned by ^(A) = Ax/Ab, the rest frame optical 
extinction can be written in terms of an observed colour, say {B — K)o, and corresponding 
intrinsic colour [B — K)i as follows: 


Av = 


e(A 


V) 


^^(A_b/1 -tr z) — ^(Ai^/1 -|- z) 


[{B - K), -{B- Kf 


( 11 ) 


where z is the source redshift. We have used the analytical fit for ^(A) as derived by Pei 
(1992) for diffuse galactic dust in the range 500A ^ A ^ 25/im. 


The rest frame extinction Ay for each source was estimated by computing the average 
of the two extinction values obtained independently from the two observed colours Bj — I 
and Bj — K^- As shown in Figure there is a tight correlation between these estimates for 
Ay, suggesting each reddening indicator may equally provide a measure of the extinction. 
From the scatter about the line of equality, the values are consistent to within ±0.5 mag. 
This correlation suggests that the B — K and B — I colours of Par kes-quasars vary with 
each other in such a way implying a characteristic spectral curvature from B to K that 
is consistent with reddening by dust. This is unlikely to arise from intrinsic relationships 
between these bandpasses in the source emission, such as a correlation between the I and 
K bands. Appropriate hne tuning would be required to reproduce the result in Figure |^. 


EDITOR: PLACE FIGURE | HERE. 
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4.2. Results and Model Comparisons 

Using the above formalism to convert the observed quantities into the source rest 
frame, Figure ^ shows the rest frame spectral indices asx and axx as a function of the 
optical extinctions. The triangles represent 2a lower limits on these indices from the 
X-ray non-detections. In the lower left corners of each figure, we show two ‘conservatively’ 
calculated error bars for both the Ay and spectral indices. From the dispersions in our 
assumed intrinsic colours (eq. [^) and scatter in Figure ^ we have typically: a(Ay) ~ 0.5 
mag. The spectral index error bar assumes a maximum error in the X-ray flux of 25% 

(ie. cr(log/ 3 ;) = 0.25; see Siebert et al. 1998) and an uncertainty in both the Bj and 
bands of 0.2 mag (see § ||). Furthermore, following Siebert et al. (1998) we also included 
an uncertainty for possible variability in the non-simultaneous X-ray and Bj (or iF„) 
measurements. This assumes a variability of a(B) = a{K) = 0.3 mag, typical for one of our 
quasars with a 20 year rest-frame timescale between the optical and X-ray measurements 
(Hook et al. 1994). Combining these uncertainties, we have a{aBx) = o'(«itx) — 0.15. 

The opticaUto-X-ray slopes (Fig. |]a) and near-infrared-to-X-ray slopes (Fig. Pr) 
appear to be somewhat anti-correlated and correlated with Ay respectively. We have 
formally computed the probabilities that the observed correlations are real by taking into 
account all lower limits on the spectral slopes and using the techniques of survival analysis. 
The correlation and regression analyses were performed using the ASURV package (Version 
1.3; La Valley et al. 1992), which is particularly designed for censored data, implementing 
the methods presented in Isobe et al. (1986). For the correlation analysis we applied the 
generalised Kendall’s tan test and for the regression analysis we used the parametric EM 
algorithm. For asx vs. Ay (Fig. |^a) we hnd that the probability for no correlation is 
P = 0.0006, i.e. the hypothesis that these quantities are uncorrelated is rejected at the 
99.94% confidence level. For axx vs. Ay (Fig. |b), we find P < 0.0001, the correlation 
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is significant at the > 99.99% confidence level. The solid lines in Figure | are our best 
regression hts to the data. Using ASURV, we have the following best regression line hts: 

aBx = (1.41 ± 0.02) - (0.08 ± 0.02) x Ay (12) 

okx = (1-18 ± 0.03) + (0.10 ± 0.02) x Ay. 

The observed anti-correlation and correlation involving obx and respectively 
appear broadly consistent with the predictions of our simple model in Figure |l|. The 
predicted ranges in these slopes as a function of the optical extinction (taking account 
uncertainties in dust-to-gas ratios and intrinsic slopes; see § |^) are represented by the 
regions within the dashed curves in Figure For Ay ~ 0, we see there is very good 
agreement between the observed slopes of Parkes quasars and our assumed range in 
intrinsic (unabsorbed) values (eq. [0) from studies of optically-selected quasars. The 
optical (and near-IR)-to-soft X-ray continuum slopes of the bluest Parkes quasars therefore 
show relatively small scatter similar to those of optically-selected quasars, strengthening 
the claim that such sources are unbiased with respect to reddening by line-of-sight dust. 
The increased scatter in these slopes when all Parkes quasars are considered can then be 
attributed to dust extinction. In particular, the range in extinctions: 0 < Ay < 6 mag, 
are consistent (within our conservative errors) with the maximum extinction ~ 4 mag 
derived by Siebert et ah (1998) from the dispersion in opticaUto-X-ray flux ratio alone. 

Our results are also consistent with studies of the optical-to-soft X-ray continua of 
a large sample of medium to moderately-hard X-ray selected AGN by Puchnarewicz et 
ah (1996). Their study suggests moderate absorption by dusty gas with approximately 
galactic dust-to-gas ratio and columns A^h < 5 x lO^^cm”^. Using the lower limit 
Nb/EB- y — 3.3 X 10^^cm“^mag“^ from equation (^) and a galactic extinction curve, this 
corresponds to ^6 mag, entirely consistent with the range found in our sample of 
radio-selected quasars. 
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EDITOR: PLACE FIGURE | HERE. 

To summarise, we have found significant correlations between optical (and near-IR)-to- 
soft X-ray continuum slope and optical extinction that are consistent with the predictions 
of a simple dust model. The results of Figure imply that if the observed scatter in optical 
slopes of Parkes-quasars is due to extinction by dust in their environs, then soft X-rays are 
absorbed at a level fully consistent with this hypothesis. 


5. Discussion 

The anti-correlation involving asx and optical extinction (Fig. |]a) can be easily 
explained by gas-dust absorption whereby optical flux is relatively more absorbed than that 
at IkeV. This anti-correlation is qualitatively similar to that claimed by McDowell et ah 
(1989) for a sample of quasars which appeared to show a weak “big blue bump” feature 
relative to their near-infrared and soft X-ray emission. These authors however interpreted 
the trend as due to an intrinsically varying blue-bump spectrum and not extinction by dust. 
A changing blue luminosity may thus mimic a variation in optical extinction. It is therefore 
possible that the observed correlation in Figure ^ja is spurious and due to a secondary 
trend, ie. from a dependence of both asx and Ay (or optical slope) on R-band luminosity. 
A direct investigation of the relationship between these quantities however shows that 
such an effect is unlikely to explain the observed trend. Also, since the optical-UV flux 
in the blue-bump component also provides a source of photoionizing flux for emission line 
regions, this model also predicts to first order a correlation between the equivalent widths of 
emission lines and optical-UV continuum slope. Our analysis of the emission line equivalent 
widths for Parkes quasars (Masci 1997, Francis et al. 1998) however strongly argues against 
this possibility. An intrinsically varying optical/UV continuum will also require the X-ray 
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continuum to vary simultaneously to satisfy the correlation involving axx in Figure |b. 
This however would be inconsistent with the trend in Figure |^a. Such a mechanism 
therefore requires physically unmotivated hne tuning of the near-infrared, optical and soft 
X-ray emission to explain the observed correlations. 

The correlation involving axx (Fig. |b) indicates that all of the optical-to-lkeV 
continuum flux appears absorbed relative to that in the near-infrared. It is encouraging to 
hnd that such a correlation is qualitatively similar to that claimed by Ledden & O’Dell 
(1983) between radio-to-soft X-ray and radio-to-optical slope for several optically reddened 
radio quasars. They concluded that absorption by associated gas and dust was the most 
likely explanation. Their statistics however were too low from which to draw any hrm 
conclusion. 

The indirect conhrmation for soft X-ray absorption from Figure ^ implies that the 
soft X-ray, optical and possibly also the near-IR continuum emission must traverse the 
same dust component. As suggested by previous studies to explain the difference in soft 
X-ray spectral properties of radio-loud and radio-quiet quasars, a signihcant fraction of the 
soft X-ray flux in radio-quasars is believed to arise from synchrotron self-Compton (SSC) 
emission associated with the radio jet emission (eg. Wilkes & Elvis, 1987, Ciliegi et ah 
1995). According to standard unihed models for AGN, the radio jet emission may extend to 
distances lOkpc from the central engine. Thus, if the soft X-rays are mostly SSC in origin, 
then the results of Figure ^ require the absorbing medium to extend on a scale similar to 
that of the radio emission. If the absorbing medium however were situated close to the 
central AGN so that the SSC X-rays suffered minimal absorption (with effectively tx ~ 0) 
then the predictions of Figure ^ will change drastically: a much steeper anti-correlation in 
Figure ^ and the opposite behaviour (ie. anti-correlation) in Figure ^3 would be expected. 
More direct studies of X-ray absorption, preferably via high quality spectral observations 


are required to further explore this issue. 


6. Conclusions 

We have tested the dust-reddening hypothesis to explain the relatively large dispersion 
in optical-to-near-infrared colours of quasars in the Parkes sample of flat-spectrum radio 
sources by searching for associated absorption by gas at soft X-ray energies. We have used 
broadband (0.1-2.4 keV) ROSAT All-Sky Survey data and pointed PSPC observations 
provided by Siebert et al. (1998), and contemporaneous optical-near-infrared photometry 
available for 119 of the 323 sources in this sample. 

A soft X-ray absorption signature was searched for indirectly by exploring the optical 
(and near-infrared)-to-X-ray continuum properties as a function of optical extinction, 
where specihc strong correlations are expected under a simple gas-dust absorption model. 

Our main conclusion is that signihcant correlations (at > 3a conhdence) are observed 
between rest frame optical (and near-infrared)-to-lkeV continuum slope and optical 
extinction that are consistent with the predictions of a simple model. This model assumes 
the range in dust-to-gas ratios and metal abnndances derived empirically from the galactic 
ISM, intrinsic (unabsorbed) slopes as observed in optically-selected quasars, and as 
snggested by previous studies, that all absorption is intrinsic to the qnasars. Under these 
assnmptions, we conclnde that soft X-rays are absorbed at a level consistent with the range 
of extinctions 0 < Ay < 6 mag observed. 

The dust associated X-ray absorption is therefore consistent with physical conditions 
(eg. gas ionization state) fonnd in the diffnse local ISM. No warm (ionised) absorption is 
reqnired, since the data does not indicate excess X-ray absorption relative to that in the 
optical compared to the galactic model predictions. This is contrary to nnmerous previons 
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studies of the soft X-ray properties of AGN. High quality X-ray spectra to detect heavy 
metal absorption edges however are required to place stronger constraints on the physical 
properties. These data will also be necessary to provide ‘direct’ estimates of the amount 
of X-ray absorption and for constraining the geometry of absorbing/emitting regions. Due 
to the faintness of many of the reddest sources, spectral data is currently unavailable, and 
requires the high signal-to-noise and resolution capabilities of future X-ray missions like 
AXAF and XMM. 

We are grateful to Siebert et ah for making available the ROSAT data (currently in 
press) and to Paul Francis for providing the unpublished contemporaneous photometry (to 
appear in Whiting et ah 1998; in preparation). 
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Fig. 1.— (a) Rest frame Rj—to-lkeV power-law continuum slope (/i. oc and (b) 

to-lkeV continuum slope as a function of optical extinction Ay predicted assuming the 
empirical ranges in galactic dust-to-gas ratio, extinction measures and intrinsic slopes as 
discussed in § Dashed lines represent mean predictions. 

Fig. 2.— Rest frame extinctions as determined from the observed Bj — I colours versus those 
determined from Bj — Kn- The diagonal line is the line of equality. 

Fig. 3.— (a) Rest frame Rj-to-lkeV continuum slope (/j^ oc and (b) Kn-to-lkeV 

continuum slope as a function of rest frame optical extinction for Parkes quasars. Triangles 
represent 2a lower limits from the X-ray non-detections. Errors in both Ay and slopes are 
represented by conservatively calculated error bars shown in the lower left, as discussed in 
§ |4.2| . Dashed lines represent the extremities of the range predicted by our simple model (see 
Fig. 13). Solid lines are best line regression hts. 











